Axion field theory and anomalous non-dissipative transport properties of 
(3+l)-dimensional Weyl semi-metals and Lorentz violating spinor electrodynamics 
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From a direct calculation of tire anomalous Hall conductivity and an effective electromagnetic 
action obtained via Fujikawa's chiral rotation technique, we conclude that an axionic field theory 
with a non-quantized coefficient describes the electromagnetic response of the (3 + f )-dimensional 
Weyl semi-metal. The coefficient is proportional to the momentum space separation of the Weyl 
nodes. Akin to the Chern-Simons field theory of quantum Hall effect, the axion field theory violates 
gauge invariance in the presence of the boundary, which is cured by the chiral anomaly of the 
surface states via the Callan-Harvey mechanism. This provides a unique solution for the radiatively 
induced CPT-odd term in the electromagnetic polarization tensor of the Lorentz violating spinor 
electrodynamics, where the source of the Lorentz violation is a constant axial four vector term for 
the Dirac fermion. A direct linear response calculation also establishes anomalous thermal Hall 
effect and a Wiedeman-Franz law, but thermal Hall conductivity does not directly follow from the 
well known formula for the gravitational chiral anomaly. 



I. INTRODUCTION 



Usually topological states possess gapped spectrum 
in the bulk. The charge and the spin quantum Hall 
states in (2-|-l)-dimensions, and the strong topologi- 
cal insulator in (3-1-1) dimensions are some of the well 
known examples of such gapped topological states of mat- 
ter. Topological considerations have also been applied 
to (3-hl)-dimensional Weyl fermions, which describe a 
gapless, semi-metallic phase with nontrivial momentum 
space topologyji In the simplest realization of Weyl semi- 
metal, there are linearly dispersing, 2-component mass- 
less Dirac fermions around a pair of diabolic points in 
the momentum space, described by the Hamiltonians 
h± = ±t;i?cr. (k—K±), where vp is the Fermi velocity, and 
(T is the vector formed out of three Pauli matrices. The 
two component Hamiltonians h± respectively describe 
the source (monopolc) and the sink (anti-monopole) of 
the Berry curvature in the momentum space. The ± 
signs also correspond to the right and the left chirality of 
the Weyl fermions. For Ak = K_|- — K_ 7^ 0, the right 
and the left Weyl fermions are separated in the momen- 
tum space and the semi-metal is topologically non-trivial. 
In contrast Ak = K+ — K_ = corresponds to topolog- 
ically trivial, 4-componcnt massless Dirac fermion. A 
Weyl semi-metal as described above, breaks time rever- 
sal symmetry and gives rise to anomalous charge and 
thermal Hall conductivities, and also possesses surface 
states with chiral dispersion in the directions orthogonal 
to Ak,^ 

Recently the (3-1-1) dimensional Weyl semi-metal 
phase has been proposed for a variety of condensed 
matter systems, such as antiferromagnetically ordered 
pyrochlore iridates^, topological insulator-normal in- 
sulator multi-layers^"—, magnetically doped topological 
insulatori^ etc. In these systems electric charge is a 
conserved quantity, and the hallmark of topology can 



be found in anomalous electrodynamic properties'"—'^"—. 
For example (3-|-l)-dimensional chiral anomaly (also 
known as Adler-Bardcen-Jackiw anomaly), which oc- 
curs for non-orthogonal arrangement of electric and 
magnetic fields may lead to a large negative magneto- 
resistance^Tii^. An anomalous charge Hall effect has also 
been predicted'"— i' when the crystallographic symmetry 
does not cause any subtle cancelation of the Berry's phase 
among different pairs of Weyl fermions^. On the other 
hand, Weyl quasiparticles can also occur for superflu- 
ids and superconductor s ^ 1 ^ 1 . For Weyl quasiparticles 
of superfluids and superconductors charge is not a con- 
served quantity. Therefore there will be no anomalous 
charge Hall effect. The topological property of Weyl su- 
perconductor will rather manifest through an anomalous 
thermal Hall effecti^. There is also tremendous current 
interest in high energy physics community in understand- 
ing how field theory anomalies affect the relativistic hy- 
drodynamics and various transport quantities'^"—. 

Conventionally the low energy, universal physics of a 
gapped topologically ordered state is described by a topo- 
logical field theory, which involves the conserved cur- 
rents and the associated gauge fields. The topological 
field theories possess quantized coefficients, and describe 
the quantized non-dissipative response of the conserved 
charge. For example the charge quantum Hall effect in 
(2+l)-dimensions is described by a Chern-Simons term 
for the electromagnetic fiel d^^i^" 



S[A] = ^J d'xe^^^A.F,, 



(1) 



where Gxy is the quantized Hall conductivity and yl^ and 
JvA are respectively the electromagnetic vector poten- 
tials and the field strengths. By varying S[A\ with re- 
spect to Ai and A2, wc obtain the quantized antisym- 
metric tensor of the Hall conductivity, and a subsequent 
variation with respect to Aq leads to rt = OxyB, which 
in turn provides the thermodynamic Streda formulai^i 
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<^xy = dn/dB, where n is the density of the matter field. 
In the presence of a boundary, the Chern-Simons action 
becomes gauge non-invariant and the offending contribu- 
tion due to the boundary is a (l-l-l)-dimensional term, 
which is proportional to the applied electric field along 
the boundary (cx axyFo2/'2'!T = axyEy/2TT, if the bound- 
ary is along y-dircction) . This gauge non-invariant con- 
tribution due to the boundary is precisely canceled by the 
chiral anomaly term of the gapless, chiral edge states via 
Callan- Harvey mechanismi ^"'^^ This is how a topological 
field theory captures (i) the quantized response of the 
conserved charge in the bulk and (ii) the bulk-boundary 
correspondence. 

For a (3+l)-dimensional, time reversal invariant, 
strong topological insulator the electromagnetic response 
is described by the axion term^^ 

S[A] = -^1 d^xef^-'P^T^^Tpx (2) 

where ~ O^ir respectively describe the quantized 
magneto-electric coefficients of a topologically trivial and 
a non-trivial insulator. In the presence of a boundary 
between the topological and the trivial insulators, the 
axion angle 9 discontinuously changes from tt to 0, and 
gives rise to a gauge non-invariant boundary term. The 
boundary term in this case becomes a Chern-Simons term 
with axy = e^/47r, and is precisely canceled by the parity 
anomaly term^iiS^ of the (2-|-l)-dimensional, two com- 
ponent Dirac fermions on the surface, via Callan-Harvcy 
mcchanism22i2^. The above action describes a quantized 
magnetoelectric effect, which manifests as a half- integer 
quantum Hall effect of the two component Dirac fermions 
on the surfaco22,. 

Therefore it is natural to ask if there is a topological 
field theory description of a Weyl semi-metal. In the lan- 
guage of relativistic field theory, the momentum space 
separation between the left and right Weyl fermions is 
captured by a purely space-like axial vector, and an axial 
vector breaks Lorentz invariance and CPT symmetry2£. 
The consequences of Lorentz violation due to CPT odd 
axial vector have been widely studied in high energy 
physics22r— . On the symmetry grounds it is well known 
that the spinor field can radiatively induce a CPT odd ax- 
ion term in the vacuum polarization tcnso r^^'^^" — . How- 
ever the precise coefficient and the expression of the axion 
angle has been a subject of dispute. 

Recently in Ref. l3lL Zyuzin & Burkov have applied 
Fujikawa's chiral rotation method^ for deriving effective 
electromagnetic action of a Weyl semi-metal and have 
found a spatially varying axion angle 9 = Ak.x, where x 
describes the three spatial coordinates in the bulk. For 
an axial 4-vector 6^ and zero Dirac mass, the induced 
axionic action is given by 

[A = ^ j e"''"^ J-^, (3) 

From this effective action they have reproduced the 
anomalous Hall conductivity a^^ = e^/(27r^)e'-'''"6i, and 



the chiral magnetic conductivity (Jch = e^/(27r^)&o- How- 
ever all the subtle aspects of chiral anomaly have not been 
addressed. In Ref. |3J, without dealing with the chiral 
anomaly, the bulk-boundary correspondence for anoma- 
lous Hall conductivity has been invoked to fix the ambi- 
guity of the axion term for a general axial 4-vector in the 
presence of a Dirac mass. However these papers have not 
answered the most important question: why do the avail- 
able calculations of vacuum polarization tenso r^^i^^" — 
fail to obtain the correct coefficient of the axionic term? 

In this paper we thoroughly investigate various as- 
pects of chiral anomaly and the axionic electrodynamics 
of (3-|-l)-dimcnsional Weyl scmi-mctal. When the exter- 
nal electric and magnetic fields are orthogonal, there is 
no (3-fl) dimensional chiral anomaly. In this case we 
show that the (1+1) dimensional chiral anomaly associ- 
ated with the surface states precisely cancel the gauge 
symmetry violating boundary term arising from the ax- 
ionic action via Callan- Harvey mechanising^, and the net 
theory indeed remains anomaly free. When the electric 
and the magnetic fields are not orthogonal or a magnetic 
field is applied in the presence of a chiral chemical po- 
tential (purely time-like axial vector), (3+l)-dimensional 
chiral anomaly shows up and this bulk chiral anomaly is 
not compensated by any surface effects. We also deal 
with the direct calculation of the CPT odd part of the 
polarization tensor and the anomalous Hall effect. From 
our detailed calculations it will become clear that, due to 
the violation of both Lorentz and spatial rotational in- 
variance by a purely space-like axial vector, any Lorentz 
or spatially rotational symmetric ultraviolet regulariza- 
tion leads to erroneous answer for anomalous Hall coef- 
ficient (inconsistent with the bulk-boundary correspon- 
dence), and renders the overall theory anomalous. Only 
when the ultraviolet regularization is in accordance with 
the reduced rotational symmetry, we obtain the correct 
anomalous Hall conductivity, and simultaneously satisfy 
the bulk-boundary correspondence. 

Our paper is organized as follows. In Sec. II we de- 
scribe the relevant Dirac Hamiltonian augmented by an 
axial vector, and discuss the emergent spectrum in detail. 
In Sec. HI wc provide the detailed solutions of the chiral 
surface states of a Weyl semi-metal and a Weyl supercon- 
ductor. We also demonstrate how (l-l-l)-dimensional chi- 
ral anomaly becomes embedded in a (2+l)-dimensional 
set up. Based on the surface state solutions we provide 
the relevant formula for anomalous charge Hall conduc- 
tivity. In Sec. IV we obtain the axion electrodynamics 
of Weyl scmi-metal by employing Fujikawa's chiral ro- 
tation technique. We show how this axionic field the- 
ory captures all the important features of the anomalous 
electrodynamics (i) anomalous Hall effect, (ii) Streda for- 
mula and (iii) bulk-boundary correspondence. In Sec. 
V we compute the anomalous Hall conductivity using 
Kubo formula, and show the detailed comparison among 
the answers obtained with different ultra-violet regular- 
izations. In this section we also show for the first time, 
that the anomalous Hall conductivity remains unchanged 
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even in the presence of a finite chemical potential, in con- 
trast to the conventional problem of anomalous Hall ef- 
fect in dilute magnetic semi-conductors^^^. In Sec VI 
we briefly address the problem of purely time-like axial 
vector or an axial chemical potential, and associated for- 
mula for the chiral magnetic conductivity. In Sec. VII 
we briefly comment on the relation between thermal Hall 
conductivity and the gravitational chiral anomaly. In 
Sec. VIII we present a summary of our results and the 
future directions, and a solution for the chiral surface 
states of a Weyl superconductor is presented in the Ap- 
pendix. 



II. MODEL HAMILTONIAN OF WEYL 
SEMI-METAL 

For simplicity we consider the following real-time 
(Minkowski space) action involving one species of 4- 
component Dirac fermion, described by 



Sn = / rfV* 



(4) 



where h and the Fermi velocity v have been set to 
unity, m is the conventional Dirac mass (Lorentz scalar) , 
and the anticommuting 7 matrices satisfy {7^,7'^} = 
2gtJ.<y^ gfiiy = is the metric tensor, 7^ = 

*7°7^7^7"^ and 5* — ^'^70. We have added a Lorentz 
symmetry violating axial 4- vector term ijih^'y^^^ip in the 
action for generality. The Hamiltonian corresponding to 
the action in Eq. 2] is described by 



H 



Sx-^^ [-i7°7^(9j + m7° + 6^7°7''7^] ^' (5) 



The complete solution for the dispersion relation in the 
presence of a general axial 4- vector and Dirac mass has 
to be obtained numerically. The exact analytical solu- 
tion can be found for (i) a general axial 4-vector and 
TO = 0, (ii) a purely space-like axial 3-vcctor and to 7^ 
(iii) a purely time-like axial vector or an axial chemi- 
cal potential and m 7^ and (iv) a light-like axial 4- 
vector and to 7^ 0. A purely space-like axial vector term 
breaks time reversal symmetry (T), but preserves inver- 
sion [V) and charge conjugation (C) symmetries, and 
consequently breaks CPT symmetry. An axial chemical 
potential 60 a-lso breaks CPT symmetry, but in a dif- 
ferent way. The axial chemical potential breaks V , but 
preserves T and C. In this section and also in the follow- 
ing Sees. III-V we will be mainly interested in a purely 
space-like axial vector, i.e, &o = 0, which is pertinent for 
the realization of Weyl semi-metal. The effects of axial 
chemical potential 69 will be discussed in Sec VI. 

In the basis of even and odd parity bands (also known 
as Dirac representation of 7 matrices). 



the Hamiltonian operator in the momentum space be- 
comes 



Hn = 



TO 1 + CT.b 

o-.k 



(T.k 
-TO 1 + CT.b 



(7) 



Thus an axial 3- vector term in Dirac basis, corresponds 
to a ferromagnetic order or uniform Zceman coupling, 
where both parity even and odd bands have identical 
g factors. If we use the chiral representation of the 7 
matrices, 









^ (tM 




^ -1 \ 




: 


, 7^' = ( 


, -a^ j 


,7^H 


, 1 



(8) 

the Hamiltonian operator in the momentum space be- 
comes 



Hch 



cr.(k + b) 

m 1 



m 1 
-cr.(k-b) 



(9) 



When m = 0, the above Hamiltonian operators de- 
scribe two Weyl fermions separated in the momentum 
space by the vector Ak = 2b, which becomes particu- 
larly transparent from Eq. [HI In this case the Hch is 
block diagonal, and chirality of the Weyl fermions after 
momentum space splitting is still determined by 7^. For 
TO 7^ 0, the eigenstates of 7^ get mixed, and the disper- 
sion relations are given by 




ik.by 



(10) 

Due to the time reversal symmetry breaking, the 
Kramers degeneracy of the conduction and the valence 
bands are removed, which is captured by the non- 
degeneracy of conduction and valence bands for a given 
value of s. The bands corresponding to s = are 
always fully gapped. But the bands corresponding to 
s = —1 are fully gapped, only when to^ > b^, and the 
system is an insulator. For to^ < b^, the conduction 
and the valence bands corresponding to s = —1, cross 
at k = zLb\/h'^ — m? , givin g rise to ri ght and left Weyl 
fermions, with Ak — 2b -^/l — m? /b"^. In the subspace 
of these emergent Weyl fermions, we can define a new 
chirality matrix 7^ , and the low energy Hamiltonian will 
acquire the form 



Tjreduced 
J^ch 



cr. ( k + b\/b^ — m? 





k — b\/b^~^rrv 
(11) 



III. CHIRAL SURFACE STATES AND 
(l+l)-DIMENSIONAL CHIRAL ANOMALY 



1 
-1 







fT-' 





,7^ = 



1 

1 



(6) 



In this section we find explicit solutions for the surface 
states of a Weyl semi-metal. To the best of our knowl- 
edge the explicit solution for the surface states of the 
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Hamiltonian in Eq.[5]has not been published in the hter- 
ature. Therefore we provide the detailed solution of the 
surface states. In addition we also demonstrate for the 
first time how (l+l)-dimensional chiral anomaly becomes 
embedded into a (2+1) dimensional setup, and provides 
a strong constraint on the bulk effective field theory. 

Consider Hch, with h = bz and a boundary in the yz 
plane. For simplicity we can model the boundary by a 
Heaviside step function of coordinate x, b{x) = b6{x). To 
obtain the surface state solution we solve the following 



It turns out that only the state corresponding to eigen- 
value +1 of the matrix 7^7^ leads to a nornializable so- 
lution. The general ansatz for such a state is 

^{x)^u{x)\^A+v{x)\-A (14) 



which is a linear combination of the two eigenstates of 
the chiral matrix 7'"'. The above ansatz for 4'{^) leads to 
the dispersion E = +ky and the functions u{x) and v{x) 
satisfy the following coupled differential equations 

[dx + k,^+ be{x)] u{x) + mv{x) = (15) 
mu{x) + [dx -k^+ be{x)] v{x) = (16) 

From these equations by eliminating v{x), we obtain the 
following second order differential equations for u(x), 

dlu{x) - {kl + m^) u{x) = 0, y X < (17) 

dlu{x) + 2b d^x) - (fcz - b^) u{x) = 0,\fx>0 

(18) 

From these equations we obtain, 

u{x <0) = A eV'^H"'^ ^ (19) 
v(x <0) ^ ~—(k, + Jkl + mA A eV'^H™' (20) 

TO V / 

u{x>{))^Be\ ) j^Ce \ ^ ' ' (21) 



v(x > 0) = 




The continuity conditions on u(x) and v{x) at x = lead 
to A = B + C and A = B + 'hZ^l^^^Q fj-Qj^ which we 
obtain C = 0, and A ~ B. Therefore the final solutions 



eigenvalue equation 

/ icr ■ W + as b 0(x) ml \ -r / \ t^t-/ \ 

/ • V7 , uof \ *(x) = -B*(x), 

(12) 

The momentum along y and z directions are good 
quantum numbers and after we substitute ^'(x) = 
^ikyy+ik^z^f^^^ in Eq. [121 wc find that the spinor ^p{x) 
satisfies 



(13) 

I 

for u{x) and v{x) are given by 

u{x < 0) = A eV^^+^ (23) 
7.(x>0)=Ae("''+^'"+"')" (24) 
v(x) = ( + Jkl + to2 ) it(a;) (25) 

For x > 0, the solution is normalizablc only when 
— — m^ < kz < — m^. These surface state solu- 
tions with E ~ +ky are responsible for producing Fermi 
arc in the ARPES measurements^. 

If we modify the Hamiltonian with a particle-hole sym- 
metry breaking higher gradient term "i'^d^'^, the en- 
ergy E can can acquire explicit dependence on k^, but 
chiral dispersion can occur only along x (in xz plane) 
and y (in yz plane) directions. Such perturbations are 
generically present in a lattice realization of Weyl semi- 
metal, and their qualitative effects in the low energy 
and momentum limit will be to modify E = +vpky to 

E w 2mF~^'^P^V' '^^^^^ have restored the Fermi veloc- 
ity of the bulk-quasiparticles, and to* is a non-universal 
effective mass. In the presence of an electric field along y 
direction, the chiral dispersion along y leads to U{1) chi- 
ral anomaly. For any allowed value of kz, wc have chiral 
anomaly of one dimensional states j^'nE^jy^ and the 
net contribution from the entire range of —^/h^ — < 
kz < Vb'^ — m^ is given by e^Ex/y/^i: J dkz/{2TT) = 
e^/(27r^)-\/6^ — m^E^jy. Thus (l-l-l)-dimensional chiral 
anomaly is embedded in a (2+l)-dimensional setup, and 
wc can express this as the following effective action 

STn^"'' = j dtdydzf{y,z,t)Ey{y,z,t) 

(26) 

where the function f{y,z,t) is related to a gauge trans- 
formation performed on the vector potential (which is 
also a chiral gauge transformation in this case). Under 
the variation with respect to /, we will obtain the anoma- 
lous Ward identity for the charge current (which is also 



iaidx - cr2ky - a^kz + a^b 9{x) 



ml 



-i CTi dx 



ml 
a2ky + a^kz 



(73 b e{x) 



V'(x) = Eipix). 
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the chiral current in this case). The presence of the chi- 
ral surface states and the associated U(l) chiral anomaly 
suggests that the bulk electromagnetic action contains 
some kind of topological term, which in the presence of 
a boundary can cancel the chiral anomaly to make the 
overall theory anomaly free^^. Akin to the edge modes of 
the quantum Hall problem, the chiral surface states give 
rise to the following anomalous Hall conductivity^!^ 

(r,y = ^Vb^-m^ (27) 



IV. AXION ELECTRODYNAMICS IN THE 
BULK VIA CHIRAL ROTATION 

In this section we describe the bulk electromagnetic 
action obtained via Fujikawa's chiral rotation technique 
in Ref. Isil For simplicity we choose m = 0, and two 
Weyl fermions are now located at k = ±b. The relevant 
fermion action coupled to electromagnetic gauge fields is 
given by 

(28) 

We can perform a chiral rotation of the spinor fields by 
an angle 9 = b.x, which removes the axial vector bilin- 
ear from the fermion sector in the action S[i!, b,A] in 
Eq. [28l Under the chiral rotation the spinor fields trans- 
form to ^ — exp(ib.X7^)vI' and ^ — J- exp(ib.x7'^). 
However, due to the non-invariance of the path integral 
measure under such a chiral rotation, the effective action 
acquires a contribution from the Jacobian of the trans- 
formation. This procedure leads to the following effective 
action 

S[^,^,b,A]^ J d''xWi-,''{d^ + iA^)^' 

+ d'x b.x e^"^"^ J-^, J-p;, (29) 

The primed spinor fields represent topologically trivial 4- 
component masslcss Dirac fermions, and the entire effects 
of the axial vector is now captured by the axionic term 
for the electromagnetic fields, 

AS^^J d^x b.x e'^'^P^ T^, J-pA (30) 

This term does not reflect chiral anomaly in three- 
dimensional bulk. Rather this term is analogous to 
the condensation energy derived for a broken symmetry 
phase. It is more useful to write the gauge covariant form 
of this action 

= ^ (1 - ^y^-^) J e'^""' A dpAx (31) 



obtained via integration by parts. The variation of A5 
with respect to the vector potential provides the differ- 
ence between the currents in the topological and the triv- 
ial sectors. Since topologically trivial 4-component mass- 
less Dirac fermion does not have any anomalous trans- 
port property, the variation of A5 uniquely specifies the 
anomalous transport properties of the Weyl fermions. 
For b ~ hz, by varying the above action with respect to 
Ai and we obtain the Hall conductivity Uxy = 
and a subsequent variation with respect to leads to 
the Streda formula^! axy = dn/dB. In the absence of a 
boundary the axionic action is gauge invariant. However 
in the presence of a boundary (say b(x) = b9{x)z), this 
action is not invariant under the gauge transformation 
Afj,{x) — >■ Afj,{x) — dfj,f{x), and there is a boundary term 

ASboundary = j 5{xi) f (x) E^lP^ OpAx (32) 

which is precisely canceled by the chiral anomaly from 
the surface states presented in Eq. [551 For to ^ 0, we 
can integrate out the higher lying bands corresponding 
to s = +1 described in Eq. [TOl and use an effective ac- 
tion corresponding to the Hamiltonian operator shown 
in Eq. 1111 Then our results for m = problem, can be 
taken over by replacing b — > Vb^ — m? . 

In high energy physics litcratur o^^i^^'^° it has been ar- 
gued that a detailed knowledge of the high energy sec- 
tor (beyond QED) or a physical input is required to 
fix the coefficient of the axionic term. However, in the 
present case, notice that the (l+l)-dimensional chiral 
anomaly of the surface states docs not depend on the 
precise ultra-violet form of the Hamiltonian. Rather 
the emergence of the chiral surface states is precisely 
connected to the low energy physics of the bulk Weyl 
fermions. Thus in the present case the condition of an 
overall anomaly free theory or the bulk-boundary cor- 
respondence, unambiguously fixes the coefficient of the 
axionic action for the case of a purely space-like axial 
vector to be 6^/(47r2) (1 - (S^^q) (see Eq. EJ)- For the 
corresponding massive theory, the coefficient becomes 
6^/(4^2) (1 _ ^l_rnV|b|20(|b|2 - to2). 



V. ANOMALOUS HALL CONDUCTIVITY VIA 
KUBO FORMULA AND ULTRA-VIOLET 
REGULARIZATIONS 

In this section we will calculate the anomalous Hall 
conductivity using Kubo formula, which requires a di- 
rect calculation of the CPT odd part of the vacuum po- 
larization tensor. In order to succinctly describe all the 
subtleties involving ultraviolet regularization, we again 
choose the simplest case of m = 0. The origin of the 
anomalous term in the electrodynamic action is rooted 
in the trace identity rr[7^7''7''7^7'^] = -Aie^^P^. The 
real-time fermion propagator corresponding to the action 
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in Eq. |4l for m = is given by 



G{k) 



{k - 5)2 



(fc + 5)2 



After substituting this propagator in the expression of 
the polarization tensor 



Tr[j''G{k)j''G{k + p)] (34) 



and evaluating the trace over the 7 matrices, we obtain 
the CPT-odd anomalous part of the polarization tensor 



(2^ 

(k - b). 



(fc + b), 



{k + bY{k + b + pY 

(35) 



(fc-6)2(fc-6+p)2 

The integral produces a term proportional to bp. It is 
important to take the momentum cut-off along the di- 
rection of the axial vector to infinity, only at the end of 
the calculations. Otherwise a naive shift of fc^ will make 
this integral to vanis h^^i^^ . 

If we now set the external momentum p in the inte- 
grand to be zero, and perform a Wick rotation to Eu- 
clidean space, and perform the fc-integral in a Lorentz 
invariant way with a 4-momentum cutoff A, we obtain 



t'^P^ bp / k^ dk I siv^eidOi 
Jo Jo 

(36) 



P^o px (27r)452 JO JO 

^^^^^'^^^X ' {k-^ + b^ + 2kb COS0i)2 



^2^ °P 



(37) 



This gives rise to an anomalous Hall conductivity and a 
chiral magnetic conductivity, which are smaller than the 
correct answers by a factor of 4 (sec below Eq. [31). 

For the calculation of Hall conductivity using Kubo 
formula, we will choose the finite temperature fcrmion 
propagator and also consider the finite density effects 
by incorporating a non-zero chemical potential /j, in the 
propagator. We further choose a pure space-like ax- 
ial vector bf^ ~ (0,b), and the external Euclidean mo- 
mentum p^ = iojn5p,fi. After performing a Matsubara 
sum and an analytic continuation to the real frequency 
iujn — > w + the dynamic Hall conductivity is found to 
be 



d^k 
(2^ 



(k + b). 



k + b| (4(k + b)2 - w2) 
b| + /i) 



tanh ^ (|k + b| - /i) + tanh ^ (|k 



+ (b ^ -b) 



(38) 



Since we are concerned with a non-dissipative conduc- 
tivity, we obtain the same answer for the zero tempera- 
ture dc Hall conductivity for both limits lu/T ^ and 
uj/T — > 00. The zero temperature dc Hall conductivity 
for /i = is given by 



d^k 
(27r)3 



(k + b); 
|k-hb|3/2 



- (b ^ -b) 



(39) 



If we perform the 3-momentum integral in a rotationally 
symmetric way, we obtain as in Ref. [2^ 



e'^'h, / k^dk smOde 



b2(27r)- JO JO 

(b2 + A:|b|cos6i) 
(fc2 + b2 + 2fc|b| cos 61)3/2 

2e^ M 
|b|3(2 



— ^e'J'bi / k^dk9{\h\-k) 
1") Jo 



6772 



(40) 



(41) 



which is three times smaller than the correct answer (see 
below Eq. [3]). We notice that both the Lorentz invari- 
ance and the spatial rotational invariance are broken by 
b. Therefore, the above two methods of regularization 
do not conform to underlying symmetry and lead to in- 
correct answers. Below we use a regularization that is 
consistent with the broken symmetries in the presence 
of b, and obtain the correct result for anomalous Hall 
conductivity using the Kubo formula. 

There is a rotational symmetry only in the plane per- 
pendicular to b, and if we perform the integral over k in 
a cylindrically symmetric way, it will be consistent with 
the underlying symmetry. We shall choose a cut-off A 
along the b direction, and keep the cut-off in the per- 
pendicular plane to be 00. With this regularization we 
obtain 



2 i-oo pA p2Tr 

j7rTE^'''^bi / k^dk^ / dfc|| / 

(27r)3 7o ^-A Jo 

fc|| + |b| 
[^i + (fc|| + |b|)2] 



d(j> X 



3/2 



(42) 



f dk\\ sgn(fc|| + |b|) 

T^r J-A 



(2 



(43) 



which is indeed the correct answer. Therefore we need 
to choose (i) a method of integration which is consis- 
tent with the reduced rotational symmetry; (ii) a finite 
cut-off along the axial vector, which can be sent to in- 
finity at the end of calculations; and (iii) we can keep 
the cut-off's in the directions perpendicular to the axial 
vector to be infinity without any trouble. In the presence 
of a non-zero Dirac mass, the expressions become more 
involved. However, this physically motivated regulariza- 
tion does yield the correct anomalous Hall conductivity 
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2V(27r2)Vb2" 



^e^^^bi for > m^. For < 



the anomalous Hall conductivity vanishes, which is con- 
sistent with the absence of the chiral surface states, and 
certainly is the correct answer for the bulk insulator. 

Now we consider the effects of finite density (n ^ 0) 
on the anomalous Hall conductivity. For concreteness we 
choose /i > 0. The zero temperature dc Hall conductivity 
is then given by 



(27r)3 
fci 



P-OO 










J -A Jo 



HI 



[ki + {^ + \h\rf' 

r.A+|b| 



bi 



X dx 



-A+|b| 



d(t) X 

(fc|| + |b|)2-^2) (44) 



(45) 



This is a remarkable result that the anomalous Hall con- 
ductivity remains unchanged at finite density. This is 
in accordance with the general principle that U(l) chi- 
ral anomaly remains unaffected by chemical potential or 
finite temperature^!^. 



VI. CHIRAL MAGNETIC CONDUCTIVITY IN 
THE PRESENCE OF AXIAL CHEMICAL 
POTENTIAL 

For an axial chemical potential 5o 7^ 0, the number of 
left and right fermions will be different, and the (3-1-1)- 
dimcnsional bulk chiral anomaly can manifest in the pres- 
ence of a magnetic field. When 60 7^ and to = 0, 
left and right Weyl fermions respectively produce parti- 
cle and hole like Fermi surfaces (which will be unstable 
against interaction effects, due to the perfect nesting of 
the Fermi surfaces) to maintain charge neutrality. For fi- 
nite chiral chemical potential, theory is truly anomalous 
and the bulk (3-|-l)-dimensional chiral anomaly in the 
presence of an electromagnetic field can not be removed 
by any surface state. It is also important to note that 
for a purely time-like axial vector and m = 0, there is no 
surface state solution. By chiral rotation method we shall 
find an axion angle 9 = bot, which leads to an equilibrium 
current along the direction of the applied magnetic field 
ji = e^5oi3i/(27r^), and the chiral magnetic conductivity 
defined by CTc/i = ji/Bi = e'^bo/i27r'^). 

For a purely time-like axial vector, we have spatial 
rotational symmetry (but no Lorcntz symmetry). Our 
method of regularization for the direct calculation of the 
chiral magnetic conductivity will be to choose: (i) a fre- 
quency cut-off A, (ii) keep the spatial cut-off to be in- 
finity, and (iii) perform the spatial momentum integral 
in a rotationally symmetric manner. This method again 
produces the correct answer for the chiral magnetic con- 
ductivity e^6o/(27r^). In the presence of a non-zero Dirac 



mass and an axial chemical potential, the dispersion re- 
lation is given by 



£^^=±1 =±v/(|k|+s&o)2 



(46) 



which describes an insulating bulk. Due to the broken 
inversion symmetry the Kramers degeneracy is absent. 
With our regularization, the chiral magnetic conductiv- 
ity for a massive Dirac fermion in the presence of an 
axial chemical potential vanishes. Therefore we find, the 
chiral magnetic conductivity is finite for massless Dirac 
fermion in the presence of an axial chemical potential, 
while for the corresponding massive theory it vanishes. 
In Refs. MM and|3l|, the chiral magnetic conductivity 
of the massive problem has been found to be e^5o/(27r^). 
However in these works, the frequency cut-off has been 
taken to be infinity and the spatial cut-off has been kept 
finite. Hence their regularization procedure is opposite 
to the one we have used. 



VII. THERMAL HALL CONDUCTIVITY AND 
GRAVITATIONAL ANOMALY 

Since an axial 3-vector breaks time reversal symme- 
try, on symmetry ground we expect an anomalous ther- 
mal Hall conductivity. A linear response calculation in 
the bulk (after accounting for diathermal contributions) 
leads to an Wicdeman-Franz law between anomalous 
charge and thermal Hall conductivitiea^i^ for a Weyl 
semi-metal. The thermal Hall conductivity k*^ is found 
to be 



TT^klT 



klT 



(47) 



This result is also consistent with the heat transported 
by the chiral surface states^i^. For a Weyl superconduc- 
tor there is no charge Hall conductivity. The anomalous 
thermal Hall conductivity for the ferromagnetic super- 
conductor of described in Ref . [13 is given by 



24 ' 



(48) 



where Akg is the momentum space separation of the 
Weyl fermions corresponding to the spin projection s. 
In contrast to the Weyl semi-metal there is an addi- 
tional pre- factor of 1/2, which reflects that we are deal- 
ing with BdG quasiparticles in the bulk^. The pre- 
factor of 1/2 is also consistent with the fact that surface 
states are now described by chiral Majorana fermions. 
Recently the gravitational chiral anomaly has been ar- 
gued to be responsible for a quantized thermal Hall ef- 
fect on the surface of a topological superconductor in 
class Dlll^-i^. However we find that the linear response 
thermal Hall conductivity is not directly captured by the 
well known gravitational anomaly formula, even though 
the gravitational anomaly formula correctly captures the 
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bulk-boundary correspondence at zero temperature, and 
the anomaly formula also verifies the correct conformal 
charge of the surface states^. Similar concerns regard- 
ing the connection between the thermal Hall effect and 
the gravitational chiral anomaly have also expressed in 
Refs.li^andlii 



VIII. SUMMARY AND FUTURE DIRECTIONS 

In this paper we have explored the relationship be- 
tween the chiral anomaly and the effective electromag- 
netic action of a (3-|-l)-dimensional Weyl semi-metal in 
detail. We have obtained explicit solutions for the chi- 
ral surface states and associated (l+l)-dimensional chi- 
ral anomaly. Based on the chiral anomaly of the surface 
states we have demanded the existence of an anoma- 
lous bulk action, such that the overall theory becomes 
anomaly free via Callan-Harvey mechanisro^S. From a 
subsequent calculation of the bulk charge Hall conduc- 
tivity via Fujikawa's chiral rotation method and also by 
Kubo formula, we have established that the anomalous 
electrodynamic properties in the bulk is described by ax- 
ion electrodynamics. Unlike the strong topological insu- 
lator, which has gapped spectrum and a quantized co- 
efficient for axionic action, the coefficient of the axionic 
electrodynamics for Weyl semi-metal is proportional to 
the momentum space separation of the Weyl nodes. The 
momentum space separation of the Weyl fermions can 
be renormalized due to disorder and interaction effects. 
However, weak disorder and long range Coulomb in- 
teractions turn out to be irrelevant perturbations^i^ 
for (3-|-l)-dimensional Dirac and Weyl fermions at zero 
chemical potential. Therefore at zero chemical poten- 
tial, the renormalization effects on the anomalous Hall 
conductivity will become unimportant at low tempera- 
tures. We have showed that a finite chemical potential 
does not change the anomalous Hall conductivity of the 
clean, noninteracting Weyl fermions. However the com- 
bined effects of disorder and interaction on the anomalous 
Hall conductivity at a finite density is an interesting open 
problem. 

We have pointed out the subtleties involved in the di- 
rect calculation of anomalous Hall conductivity by using 
Kubo formula, and also provided a rcgularization scheme 
that respects the underlying symmetry. Our rcgular- 
ization does produce the correct answer for the anoma- 
lous Hall conductivity for a Lorentz violating spinor elec- 
trodynamics, where the Lorentz violation is caused by 
a purely space-like axial vector bilinear for the Dirac 
fermions. Our rcgularization provides the correct chi- 
ral magnetic conductivity for a masslcss Dirac fermion 
in the presence of an axial chemical potential or a purely 
time-like axial vector. Based on our rcgularization we 
have found zero chiral magnetic conductivity for a mas- 
sive Dirac fermion in the presence of an axial chemical 
potential, which is at odds with the result reported in the 
litcraturei^i^'^. Finally, a linear response calculation of 



the thermal Hall conductivity has shown the existence of 
a Wiedemann- Franz law between the charge and the ther- 
mal Hall conductivities. In contrast to the recent claims 
in the literature, we find that the thermal conductivity 
does not directly follow from the well known gravitational 
anomaly formula. The detailed connection between the 
thermal Hall effect and gravitational anomaly will be dis- 
cussed in a separate publication^. 
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Appendix A: Surface states of Weyl superconductor 

For a Weyl superconductor described in Ref. [13, we 
can use the Hamiltonian operators of Sec. II, as long 
as we identify ^(x) to be a Nambu-spinor. However to 
make a direct connection we choose the ferromagnetic 
triplet superconductor described in Ref. [13, which has 
the same pairing symmetry as He'^ — A2 phase. Due 
to the ferromagnetic moment there are spin-split Fermi 
surfaces in the normal state. In addition the same spin 
pairing is chosen with the symmetry — iky, but with 
distinct amplitudes for two spin projections. The reduced 
BCS Hamiltonian is described by 

2m 



n 



BCS 



E 



= ±.k 



X c! , 
k,s — k,.s 



,Ck,, 



, {kx 



iky) 



h.c. 



(Al) 



where /ig. fci?,s = \/2mfis, and As are respectively the 
chemical potential, the Fermi surface radius, and the 
pairing amplitude for spin projection s. For each spin 
projection s. we have a pair of Weyl fermions located at 
kg = ±kp gZ. In order to obtain a qualitative description 
of the surface states we choose a spatially varying chem- 
ical potential ^s(x) = /i^ tanh(a;/^)^. Again we retain 
ky and k^ as good quantum numbers, and replace kx by 
—idx- For small energy and momentum, we can ignore 
and ky in comparison to dx and ky, but retain the 
only kz dependence through the quadratic term k^/2m 
to obtain 



2m 



tanh ■ 



T3 + {-iTldx + Tyky) 



A.. 

kF,s 

= EstI>s{x) 



Ipsix) 

(A2) 



In order to obtain a normalizable solution we need to 
choose Ipsix) to be an eigcnstate of Ty with eigenvalue 
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+1, which leads to the ansatz ip'^{x) = Us{x){l, —i), and 
the dispersion Eg = —Agky/kps- The function Ug{x) 
satisfies 

d:cUs{x) - - fis tanh ^) Us{x) = 0, (A3) 



which has the fohowing asymptotic property 



lim Us[x) ^ uo^sGxp 



A, 



2m 



X - tis\x\ 



(A5) 



and the solution is described by 



Us[x) = uo^sexp 



A, 2m 



cosh ■ 



(A4) 



The above solution Us{x) is normalizable only when 
— kp.s < kz < kp^s, and describes a Majorana fermion 
with chiral dispersion along —y direction. 
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